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ABSTRACT 

We demonstrate correlated rapid variability between the optical/UV and X-ray emis- 
sion for the first time in a soft X-ray transient, GRO J 1655-40: HST light curves show 
similar features to those seen by RXTE, but with mean delay of 10-20 s. We interpret 
the correlations as due to reprocessing of X-rays into optical and UV emission, with 
a delay due to finite light travel time, and thus perform echo mapping of the system. 
The time-delay distribution has a mean of 14. 6± 1.4s and dispersion (i.e. the standard 
deviation of the distribution) of 10.5 ± 1.9 s at binary phase 0.4. Hence we identify the 
reprocessing region as the accretion disc rather than the mass donor star. 

Key words: accretion, accretion discs - binaries: close - stars: individual: Nova Sco 
1994 (GRO J1655-40) - ultraviolet: stars - X-rays: stars 
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Soft X-ray transients (SXTs), a lso referred to as X-ray no- 
vae, (Tanaka & Shibazaki 1996) are low-mass X-ray bina- 
ries (LMXBs) in which long periods of quiescence, typically 
decades, are punctuated by very dramatic X-ray and opti- 
cal outbursts, often accompanied by radio activity as well. 
It is commonly held that the optical emission seen during 
outburst arises from reprocessing of X-rays by the accretion 
disc and/or secondary star (see King & Ritter 1998 and ref- 
erences therein). It is then natural to look for correlated X- 
ray/optical variability with a view to performing echo map- 
ping of the system, a technique which has had great success 
in the study of active galactic nuclei (Peterson 1992 ) . Such 



correlated variability has been seen with low time resolution 
in the persistent LMXB Sco X-l (Ilovaisky et al. 1980, Petro 
et al. 1981) and reprocessed optical pulsations in He r X-l 



have been used to estim ate the system parameters (Mid- 
dleditch & Nelson 1976), but echo mapping has yet to be 



fully applied to an X-ray binary. 

The SXT GROJ1655-40 was discovered in 1994 July 
when the Burst and Transient Source Experiment (BATSE) 
on GRO observed it in outb urst at a level of 1. 1 Crab in 
the 20-200 keV energy band (iHarmon et al. 19951). After a 



period of apparent quiescence from late 1995 to early 1996, 
GRO J 1655-40 went in to outburst again in late 1996 April 
( Remillard et al. 1996 ), and remained active until 1997 Au- 
gust. During the early stages of this outburst we carried 
out a series of simultaneous HST and RXTE visits. One of 
the primary goals of this project was to search for correlated 



variability in the two wavebands. The long-term evolution of 
the outburst argued against significant reprocessing, as the 
s eemingly anticorr elated optical and X-ray fluxes observed 
(Hynes et al. 1998) are not to be expected if the optical flux 
is reprocessed X-rays. Nonetheless, significant short term 
correlations were detected. In this paper we analyse these 
and use them to constrain the system geometry during out- 
burst. We will perform a more comprehensive analysis of 
other aspects of the variability observed during the outburst 
in a subsequent paper. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 HST 

Our HST observations took place between 1996 May 14 
and July 22 usin g the Faint Object Spectrograph (FOS) 
in RAPID mode ( Keyes et al. 1995| ). The full log is pre- 
sented in Hynes et al. (1998). In this work we focus on the 
most promising set of observations taken on June 8 with the 
PRISM and either blue (PRISM/BL) or red (PRISM/RD) 
sensitive detectors. This configuration delivered a series of 
time resolved spectra, termed groups, (useful spectral range 
~ 2000 - 9000 A) from which we constructed 2-3 s light 
curves. Over most of the spectral range of the PRISMS, the 
spectral resolution is very poor, hence we can only study 
continuum variations. 

We extracted the light curves from the count rates 
provided by STScI, performing background subtraction by 
hand. This is necessary as the standard background model 
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is kno wn to underestima te the backgrounds by up to 30 per 
cent (Keyes et al. 1995). We rescaled the standard back- 



ground models to match unexposed regions on a group-by- 
group basis and subtracted this to isolate source counts. 
We then integrated the counts over the source spectrum to 
obtain light curves. The effective bandpass can be defined 
by the FWHM of the count rate spectra: 3100-4800 A for 
PRISM/BL and 3800-7400 A for PRISM/RD. 

A final subtlety involves the start times of the groups. 
As discussed by Christensen et al. (|1997|), FOS RAPID mode 



may produce groups unevenly spaced in time: the "too rapid 
RAPID" problem. Because of this the standard data prod- 
ucts contain an uncertainty in the start times of individual 
groups of -0.255 s/+0.125 s. We therefore had our start times 
recalculated using the rapid_times program at STScI which 
reduces the relative uncertainty between groups to < 1 [is. 
There is still an unavoidable 0.255s uncertainty in the abso- 
lute start time of each exposure (i.e. the zero point of each 
series of groups.) 



2.2 RXTE data 

The RXTE /PCA data included in this paper were obtained 
on 1996 June and contains 4 segments of exposure about 
3.5 ks each. The original data were taken with two stan- 
dard EDS modes and three additional modes for high time- 
resolution (up to 62 fjs) studies. The standard-1 mode has a 
0.0125 s time resolution and a single band covering the en- 
tire 256 energy channels. Since the highest time resolution 
from the HST data is only about 2-3 s, the light curves for 
the correlation study performed in this paper are extracted 
from the standard-1 mode data in 1 s time bins using the 
saextrct task in the FTOOLS software package. The average 
count rate of GRO J 1655-40 during the above mentioned 
observing period is more than 2 x 10 4 counts s _1 of "good 
events" (i.e., after 80% of internal background events are 
rejected by the anticoincidence logic), thus no background 
(< 35 counts s _1 ) subtraction procedure is necessary. 

The relative timing accuracy of the RXTE data is lim- 
ited only by the stability of the spacecraft clock which is 
good to about 1 /xs or less. The absolute timing accuracy, 
however, is also limited by uncertainties in the ground clock 
at the White Sands station and other complications. For 
data taken before April 29, 1997, the absolute timing ac- 
curacy is estimated to be about 8 fjs which is substantially 
better than that of HST /FOS and certainly sufficient for 
our correlation study. 

The RXTE light curves of GRO J1655-40 exhibit vari- 
ability on various time scales. On the 1 s or less time scale, 
there are flickerings with RMS of 10% or so. On longer time 
scales up to a few hundred seconds, however, the amplitude 
of broad peaks and troughs (they are not necessarily follow- 
ing each other) can be as high as 50% which are the source 
for the induced optical variability we detect. 



3 COMPARISON OF LIGHT CURVES 

In Fig. |l] we show the light curves obtained in the June 8 
visit. Of the four visits on which both HST and RXTE 
observed GRO J1655-40, this achieved the best coordina- 
tion. Fig. ^a) shows the full data set to indicate the extent 



of simultaneous coverage. In Fig. |l|b) we show the correla- 
tions present in the third pair of light curves. To illustrate 
wavelength dependence we show both UV (2000-4000 A) 
and blue (4000-6000 A) light curves. While the main fea- 
ture around 1200 s is present in both, the smaller correlated 
features are more prominent, or only present at all, in the 
UV light curve, suggesting that the source of variability rep- 
resents a larger fraction of the total light in the ultraviolet 
than in the optical. There are also some features which are 
strong in the X-ray light curve, e.g. at 1700 s, but which do 
not appear in either HST light curve. These conclusions are 
borne out by a similar close comparison of the other light 
curves, and it is clear that the relation between X-ray and 
optical emission is complex. It may be, for example, that the 
observed X-ray variations originate from different locations, 
some of which can illuminate the disc and some of which 
cannot. 

Another example of uncorrelated variability may the 
apparent downward step in the second PRISM/RD light 
curve in Fig. [j]a). There may also be a step in the X-ray 
lightcurve, but the two do not match well, and the opti- 
cal light curve could not be reproduced as a convolution of 
the X-ray light curve with a Gaussian transfer function (see 
Sect. 4.2). The pronounced step in the optical light curve 



also results in a strong autocorrelation, which leads to a ver y 
broad peak in the cross-correlation function (see Sect. |4.l[ ). 
We therefore truncated this light curve just before this point 
to simplify the analysis in the subsequent sections. 



4 ANALYSIS 

4.1 Cross correlations 

We begin by performing a cross-correlation analysis. This 
will identify correlations and reveal the mean lag between X- 
ray and optical variability. The technique is commonly used 
in the study of correlated variability from active galactic nu- 
clei (AGN) where two methodologies have b een developed 
the Interpol ation Correlation Function, ICF ( |Gaskcll fc Pe- 



t erson 19871), and the Di screte Correlation Fun ction, DCF 



( |Edelson fc Krolik 198^ ). White & Peterson (^99J con- 



trasted the relative merits of the two and suggested some im- 
provements. We have tested both methods on our data and 
found no significant differences, so choose to adopt the ICF 
method with one important modification. This is that since 
the RXTE data has a finer time resolution than the HST 
data, we approximately integrate the RXTE light curve over 
each lagged HST timebin, rather than simply interpolating. 

We show a section of the resulting ICFs from the June 8 
visit, centred at zero lag, in Fig. ^ All four show roughly 
coincident peaks, significant at the 3cr level; if the data is 
combined to yield a single ICF, then the significance of the 
combined peak is 8a. The ICF peaks occur at lags in the 
range 12-24 s. While other peaks are seen in individual ICFs, 
since they are not repeated in more than one pair of light 
curves, we do not judge them significant. 



4.2 Fitting transfer functions 

To characterise the distribution of time delays present be- 
tween the RXTE and HST lightcurves, we fitted param- 
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Figure 1. a) HST and RXTE light curves from 1996 June 8. 
Note that the first two HST light curves use PRISM/RD, while 
the latter two use PRISM/BL. For clarity, the light curves have 
been rebinned to a time-resolution of ~10s. The relative count 
rates of RXTE vs PRISM/RD vs PRISM/BL have been rescaled 
but the zero point is correct. Typical errors are 0.7 per cent for 
the RXTE light curves, 1.5 per cent for PRISM/RD and 3.4 per 
cent for PRISM/BL. b) A close up of the third pair of light curves 
from a) showing correlated variability and colour dependence in 
the HST light curves. The zero point is arbitrary, but the vertical 
axis indicates the amplitude of variations relative to the mean 
count rate. The wavelengths compared are 2000-4000 A (UV) and 
4000-6000 A (blue). 



eterised model transfer functions. In this modelling we 
predict the HST light curve by convolving the observed 
RXTE lightcurve with a Gaussian transfer function (i.e. 
the time-delay distribution). We then judge the "badness- 
of-fit" by calculating the \ 2 over the data points in the HST 
lightcurve. 

This model adopts the measured RXTE lightcurve ver- 
batim, thus ignoring the statistical errors in the RXTE mea- 
surements. This is an acceptable approximation because the 
signal-to-noise ratio for detecting variations is much higher 
for the RXTE data than for the HST data. 

The Gaussian transfer function 



\/27tAt 



has 3 parameters: the mean time delay to, the dispersion 
or root-mean-square time delay (hereafter RMS delay) At, 
which is measure of the width of the Gaussian, and the 
strength of the response, V?, which is the area under the 
Gaussian. 

Figure bl shows the synthetic light curves from the Gaus- 



0.2 
0.1 
0.0 
-0.1 

0.2 
0.1 
0.0 

-0.1 




0.2 
0.1 
0.0 
-0.1 



Exp. 4 (PRISM/RD) 




Exp. 5 (PRISM/RD) 



0.2 \- Exp. 6 (PRISM/BL) 
0.1 
0.0 
-0.1 



Exp. 7 (PRISM/BL) 




-100 

Lag (s) 



Figure 2. ICFs for June 8 data. Dashed lines show 3<r confidence 
limits expected for uncorrelated variability. All four plots show 
features at ~ 20 s which are significant at the 3a level, although 
only Exp. 6 can be considered more than marginal at this level. 



sian superimposed over the 4 HST lightcurves. The principal 
features of the HST lightcurves are reproduced well in the 
synthetic lightcurves. 

Figure H shows the results of fitting the Gaussian trans- 
fer functions to Exp. 6. Panel (a) shows the constraints im- 
posed by the data on the mean and RMS delay. The best 
fit for the Gaussian fitting has Xmin/779 = 1.229. Here the 
2-parameter 1-sigma confidence region is defined by the con- 
tour x 2 = Xmin + 2.3, and the greyscale indicates relative 
probability. Panel (b) uses Monte-Carlo error propagation 
to indicate the range of uncertainty in the delay distribu- 
tion. This plot shows 10 Gaussians selected at random with 
the probability distribution indicated in panel (a). 

The Gaussian model investigated above is acausal be- 
cause it permits HST response preceding the X-ray driving. 
We can construct a causal model by truncating the negative 
delay tail of the Gaussian. Panels (c) and (d) of Fig. ^ il- 
lustrate the results for the causal Gaussian fits in the same 
format as panels (a) and (b) of the same figure. The best fit 
for the causal Gaussian fitting has Xmin/779 = 1.230. The 
parameters to and At become anti-correlated because the 
data to first order constrain the first and second moments 
of the delay distribution. 

Table |l| summarises the results of fitting Gaussian trans- 
fer functions to the 4 data segments. For all 4 data segments 
the mean and RMS delays are roughly consistent within 
the 1-parameter 1-sigma uncertainties. A weighted average 
yields to" = 14.6 ± 1.4 s, AT = 10.5 ± 1.9 s. The total re- 
sponse, is normalised using the mean count rate for the 
individual light curve. This normalised total response ap- 
pears to be roughly constant for the two lightcurves from 
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Figure 4. Top, acausal (left) and causal (right) x 2 surfaces for Gaussian transfer function fitting to the HST lightcurve 3. Bottom, 
best-fit (solid line) and trial (dashed line) time-delay transfer functions to HST Exp. 6. The response is iP(t). 



Table 1. Summary of results from parameterised Gaussian fitting 
to HST lightcurves. The rows give respectively the photometric 
phase, number of data points, reduced x 2 , mean delay, RMS delay 
and integrated response for each light curve. Uncertainties are 
based on 1-parameter, 1-sigma confidence regions. The integrated 
response, has been normalised by dividing by the mean of each 
HST light curve. 
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PRISM/RD, at a value of * ~ 30 x 1CT 3 and for the two 
lightcurves from PRISM/BL, at a value of * ~ 50 x 10~ 3 . 
This difference between the values of $ shows that the vari- 
ability is stronger at short wavelengths, and suggests that 
the reprocessing therefore occurs in a relatively hot part of 
the system. 



5 DISCUSSION 

It is clear from our analysis that reprocessing with a mean 
time delay of under 25 s dominates. This is the size of lag 



to be expected from the accretion disc assuming established 
system parameters (Orosz & Bailyn 1997, van der Hooft 
et al. 1998). Together with the narrowness of the trans- 
fer function (RMS delay ~ 10 s), this means that disc 
reprocessing must dominate over the secondary star, from 
which lags of greater than 40 s are expected at this binary 
phase (<f) ~ 0.40). Even allowing for the maximum uncer- 
tain ty in system parameters estimated by van der Hooft et 
al. (1998), our results are definitely inconsistent with the 
dominant source of reprocessing being the companion star. 
The response for the Gaussian fits, when normalised to the 
count rate for the light curves are consistent for the two 
PRISM/BL light curves and for the two PRISM/RD ones. 
The higher response, and hence reprocessing fraction, for the 
PRISM/BL agrees with this analysis of the mean and RMS 
delays, showing that the accretion disc is the most important 
region for reprocessing of X-rays in GRO J1655-40. 

Why is this the c ase? D uring the activity observed by 



van der Hooft et al. ( 1997 ), light curve analysis revealed 
that X-ray heating of the companion star was important; 
if that was the case here, we should see echoes originating 
from the companion star. One explanation for why we do 
not is that the X-ray absorbing material in the disc may 
have a significant scale height above the mid-plane so that 
the companion would effectively be shielded from direct X- 
ray illumination, thus reducing the strength of reprocessing. 
For this to be the case, then the shielding material must rise 
to H/R~ 0.25. 

If the accretion disc in GRO J1655-40 is significantly 
irradiated, we now must ask why the long term optical light 
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Figure 3. Best-fit synthetic lightcurves for all 4 data sets from 
Gaussian fitting, superimposed over actual HST lightcurves. The 
time axis is relative to the start of each exposure. The lightcurves 
have been rebinned by a factor of two for clarity. 



curve appears almost anticorrelated with X-ray behaviour? 
One possibility is that the irradiating X-ray flux is being in- 
creasingly attenuated by a disc corona. As the optical depth 
of the Compton scattering corona increases, we might nat- 
urally expect two consequences: firstly the hard X-rays, be- 
lieved to be produced by Comptonisation, would increase. 
Secondly, the irradiation of the outer disc would decrease, as 



X-rays moving nearly parallel to the disc must penetrate a 
large optical depth of scattering material. We would thus ex- 
pect that the optical/UV and hard X-ray light curves might 
appear anticorrelated, exactly as obser ved. Indeed, a related 
scenario was proposed by Mineshige (1994) to explain why 
the onset of an optical reflare in X-ray Nova Muscae 1991 
occured simultaneous to a decrease in the Thomson optical 
depth deduced from X-ray spectra. 

These observations represent a step forward in our un- 
derstanding of SXTs. This is the first time that correlated 
optical-X-ray variability has been detected from a transient 
source and supports the picture, often assumed, that SXTs 
in outburst have irradiated accretion discs similar to those 
in persistent LMXBs. 
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